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ABSTRACT 

We present two-dimensional gas-kinematic maps of the central region in Centaurus A. The adaptive 
optics (AO) assisted SINFONI data from the VLT have a resolution of 0'.'12 in K-band. The ionized gas 
species (Br7, [Fen], [Sivi]) show a rotational pattern that is increasingly overlaid by non-rotational 
motion for higher excitation lines in direction of Cen A's radio jet. The emission lines of molecular 
hydrogen (H2) show regular rotation and no distortion due to the jet. The molecular gas seems to be 
well settled in the gravitational potential of the stars and the central supermassive black hole and we 
thus use it as a tracer to model the mass in the central ±l'/5. These are the first AO integral-field 
observations on the nucleus of Cen A, enabling us to study the regularity of the rotation around the 
black hole, well inside the radius of influence, and to determine the inclination angle of the gas disk 
in a robust way. The gas kinematics are best modeled through a tilted-ring model that describes 
the warped gas disk; its mean inclination angle is ~ 34° and the mean position angle of the major 
axis is ~ 155°. The best-fit black hole mass is Mbh = (4.5^'q) x 10 7 M Q (3cr error), based on a 
"kinematically hot" disk model where the velocity dispersion is included through the Jeans equation. 
This black hole mass estimate is somewhat lower than, but consistent with the mass values previously 
derived from ionized gas kinematics. It is also consistent with the stellar dynamical measurement 
from the same AO observations, which we present in a separate paper. It brings Cen A in agreement 
with the Mbh — c relation. 

Subject headings: galaxies: kinematics and dynamics - galaxies:structure - galaxies: individ- 
ual(NGC 5128) - integral-field spectroscopy 
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1. INTRODUCTION 

During the last few years it has been realised that 
most, if not all, nearby luminous galaxies host a su- 
permassive black hole (BH) in t heir nuclei with masses 
in the 10 6 - 10 10 M© range (e.g. iFerrarese fc Fordll2005l . 
and references therein). The black hole mass (Mbh) is 
tightly related with mass or luminosity of the host stel- 
lar spheroid, bulge, (e.g. iKormendy fc R ichstone 119951 : 
IMarconi fc Huntll2003t lHaring fc Rixll2004D and~with the 
stellar velocity disper sion, a, (jFerrarese fc Merritt] 120001 : 



iGebhardt et al.lfeOOOf ) . These correlations have an amaz- 
ingly low scatter, perhaps surprisingly low, since the 
quantities Mbh and Mbuigc/c probe very different scales. 
These facts indicate that the formation of a massive BH 
is an essential ingredient in the process of galaxy forma- 
tion. 

The mass of the black hole at the center of NGC 5128 
(Centaurus A), the most nearby elliptical galaxy, is 
still under debate. Centaurus A hosts a powerful ra- 
dio source and an AGN revealed by the presenc e of 
a powerful radio and X-ray jet (see llsraell I1998L for 
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a rev i ew; see also iTingay et al.l 119981: Hardcastle et al.l 
120031 : iGrandi et all 120031 : lEvans et alfEoM and ref- 
erences therein). Rece nt stellar d y namic al measure- 
ments and modeling by ISilge et ahl (|2005l ) result in a 
black hole mass of 2.4xlO 8 M0 (for an edge-on model), 
while different gas-dynamical studies found masses in 
the range of 6.1 x 10 7 M Q to 2 x 1O 8 M depending 
mainly on the inclination angle of the modeled gas disk 
(IMarconi et al.ll200lLj2006HHaring-Neumaver et all2006t 
iKrainovic et al.l l2007f) . Gi ven a velocity di spersion for 
NGC 5128 of 138 kms" 1 (ISilge et al.ll2005h . we would 
expect a BH mass around 3 x 10 7 M Q from the Mbh — cr 
relation. If this BH mass is correct, NGC 5128 has 
the largest offset from the M bh ~ g relation ever mea- 
sured (taken the m ass values of lMarconi et al.l (|2001h and 
ISilge et ail (|2005h ) . Based on its stellar velocity disper- 
sion and the black hole mass expected by the Mbh — 0" 
relation the radius of influence of Cen A's black hole 
is ~ 0'.'3. Thus, ground based, seeing limited observa- 
tions are not suitable to resolve this radius. Using high- 
spatial resolution (~ O'.'l) ionized gas kinematics from 
adaptive optics assisted and spac e based observations , 
iHaring-Neumaver et al.l (|2006f ) and lMarconi et all (|2006l ) 
find masses reduced by a factor 3-4 compared to previous 
measurements. They were, however, limited to long-slit 
data at a few position angles and not able to precisely 
constrain the inclination angle of the modeled gas disk. 
With the availability of integral-field spectroscopy (IFS) 
in the near infrared (IR) the gas as well as the stars 
can be mapped in two dimensions even in dust-shrouded 
galaxy centers, like Cen A. Black hole masses can be 
derived separately from stars and gas from the same 
complete data set, and the gas geometry is constrained 
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by two-dimensional data. The star and gas results can 
then be compared to assess the reliability of the mo d- 
eling techniques. S INFONI (jEisenhauer et all l2003a Hbl; 
iBonnet et aLll2004f) at the Very Large Telescope (VLT) 
combines IFS with the resolving power of adaptive optics 
assisted observations and provides data at a spatial reso- 
lution of ~ O'.'l in K-band. For Cen A the radius of influ- 
ence of the black hole at the center should be comfortably 
resolved. Still, this resolution is far from anywhere near 
the Schwarzschild radius of the black hole and cannot 
constrain the central object to be a black hole. However, 
Cen A ranks among the best c ases for a supermass ive 
black hole in galactic nuclei (see Marconi et al. 2006 for 
a discussion.) 

For our dynamical model we assume a distance to 
NGC 5128 of 3.5 Mpc to be consistent with all previous 
mass determinations. Recent dis tance measurements are 
in the range 3.4 Mpc t o 4.2 Mpc (llsrael|1998t iTonrv et all 
120011: iReikubal 12004 iFerrarese et aiTl2007f ). with typi- 
cal uncertainties of ~ 8%. At the assumed distance of 
3.5 Mpc, 1 arcsecond corresponds to ^17pc. 

This paper is a follow-on t o the work presented in 
lHaring-Neumaver et all (|2006l ) (hereafter HN+06), us- 
ing high spatial resolution Naos-Conica (NaCo) imaging 
and spectroscopy data to get an accurate measurement 
of Cen A's black hole mass. While the analysis of HN+06 
was restricted to [Fe n] kinematics along 4 long slit po- 
sitions, we study in detail the kinematics of different gas 
species at the center of NGC 5128 in two dimensions. 
Unsurprisingly, different gas species exhibit different be- 
haviors. While the (highly) ionized gas shows (strong) 
influence by the jet (both in surface brightness and in 
the kinematic maps), the molecular gas (H2) seems to 
"feel" only gravity. This is the reason why we focus on 
H2 when we construct a dynamical model to measure the 
mass of the central supermassive black hole. Moreover, 
the two-dimensional data allow us to determine the in- 
clination angle of the gas disk in a robust way. Taken 
together, these advancements over the work of HN+06 
and other previous studies, greatly decrease the overall 
uncertainty in Cen A's black hole mass measurement. 
SINFONI stellar kinematics and black hole mass mod- 
eling are consistent with the gas-dynamical study and 
black hole mass measurement presented here, and will 
be presented in a separate paper (Cappellari et al., in 
prep.). 

The paper is organized as follows: Section 2 describes 
the observations and the data reduction. Section 3 
presents the gas morphology and kinematics and Sec- 
tion 4 our dynamical model. Section 5 gives the modeling 
results and Section 6 discusses them. 

2. OBSERVATIONS AND DATA REDUCTION 

All observations presented here were taken with SIN- 
FONI on the UT4 (Yepun) of the VLT of the Euro- 
pean Southern Observatory (ESO) at Cerro Paranal, 
Chile, on March 23 and April 1, 2005. SINFONI con- 
sists of a cryo genic near-infrared integ ral field spectrom- 
eter SPIFFI ijEisenhauer et all l2003a ilbl) coupled to the 
visible curvature ad aptive optics (AO) system MACAO 
(jBonnet et al.l 120031 ). For the Cen A nucleus, the SIN- 
FONI AO module was able to correct on an R^14mag 
star 36" South- West of the nucleus in excellent seeing of 
0'.'5, reaching nearly the diffraction limit of the telescope 



in the K band. With the appropriate pixel scale selected 
(0'.'05 x O'.'l per spatial element), the spectrograph was 
able to obtain spectra across the entire K band (approxi- 
mately A1.93-2.47/im) at a spectral resolution of R~4000 
and covering a 3" x 3" field of view, in a single shot. A 
total of 10 sky (S) and 15 on-source (O) exposures of 
900 s each followed a sequence OSO... and were dithered 
by only a few pixels to allow removal of bad pixels and 
cosmic rays. The frames were combined to make the fi- 
nal K band data cube with a total on-source exposure of 
13500 s. 

In the same way, we obtained H-band data (A1.43- 
1.87/im) with a slightly lower spectral resolution of 
R~3000 covering the central 3" x 3" with a scale of 
0705 x O'.'l per pixel. The overall integration time for 
H-band was 3600 s. 

The data were reduced using the SINFONI data reduc- 
tion pipeline provided by ESO. The nearest sky exposure 
was subtracted from the object frames, after which the 
data was flatfielded and corrected for bad pixels. Dis- 
tortion correction was done based on the position of OH 
line edges in sky frames. The wavelength calibration was 
based on the nightsky OH lines in H-band and on OH 
lines and arc lamp frames in K-band. Using the slit- 
let positions derived from the sky frames and the ob- 
tained wavelength calibration a three-dimensional data 
cube were created from each object frame. Finally the 
slitlet positions were adjusted if necessary based on cross- 
correlation with NaCo broadband images (taken from 
HN+06). The typical shift between the slitlets is in the 
range of ~ 0.1 — 0.3 pixels. This correction is applied to 
retain the full spatial resolution of the data. After cor- 
recting for telluric features and flux calibration the data 
cubes were mosaicked together. The stars Hip079775 
(B3III) and Hip68359 (B8V) were used both as telluric 
standards and flux calibrators. The flux calibration has 
an accuracy of about 5-10%. 

2.1. Spatial resolution 

To assess the spatial resolution of the data, the point- 
spread function (PSF) is estimated from the unresolved 
active galactic nucleus (AGN), as discussed in more de- 
tail by HN+06. No additional PSF calibration frames 
using stars were taken, since the PSF in adaptive op- 
tics observations depends on several factors such as at- 
mospheric conditions, brightness of and distance to the 
guide star, contrast of the guide star to the background, 
and these conditions are difficult to duplicate for the sci- 
ence and calibration observations. The most accurate 
measure of the science frame PSF in the vicinity of the 
AGN is therefore achieved directly on the unresolved nu- 
cleus. 

We describe the normalized PSF empirically by a sum 
of two Gaussian components; one narrow component de- 
scribing the corrected PSF core (o~ c ) and one broader 
component (a s ) which can be attributed to the seeing 
halo plus extended emission: 

PSF(r) = ^e-^ + ^e-^, (1) 

where F is the ratio of the flux of the narrow component 
and the total flux of the PSF (F = flux c /flux tot ai)- The 
quantity F provides a rough approximation of the Strehl 
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Fig. 1— Double Gaussian fit of the K-band SINFONI PSF 
on the continuum just next to the H2 line. The FWHM of the 
narrow and broad component is 0'.'12 and ( X'30 , respectively. The 
estimated Strehl ratio is 17% (see Section 12.11 for a definition of 
these quantities) 

ratio (S) which gives the quality of an optical system; S is 
defined as the observed peak flux divided by the theoret- 
ically expected peak flux of the Airy disk for the optical 
system (S = peakflux c /peakfluxAi r y)- For the following 
analysis it is sufficient to measure the quantity F, which 
also gives an estimate of the quality of the adaptive op- 
tics correction. 

For the K-band cube, a PSF model of the above form 
is fitted to the peak in surface brightness in the collapsed 
wavelength range of A2.00 to 2.10 /xm (just next to the H2 
line). The two components of the PSF model shown in 
comparison to the data in Fig. Q] have a full width at half 
maximum (FWHM) of 0"12 and 0'.'30, for the narrow and 
broad component respectively. The estimated Strehl ra- 
tio is 17%. The width of the broad component compares 
well to the FWHM of the seeing disk as measured by the 
seeing monitor during the observations (FWHMy ~ 0'.'5 
transformed to K-band as in HN+06, FWHM^- ~ 0'.'38). 

2.2. Subtraction of the stellar and non-stellar continua 

The total or averaged spectrum of the central 0'.'8 x 
0'.'8 of Cen A (Fig. ^ shows strong CO absorption 
lines at 2.3-2.4 fim indicating the stellar continuum. 
We use the penalized p ixel fi tting method (pPXF) of 
iCappellari fc Emselleml (|2004h to fit the stellar contin- 
uum with a positive linear combination of stellar tem- 
plates. As template stars serve six late type stars that 
were observed with SINFONI in the same setup that we 
used for the nucleus of Cen A. This ensures that the in- 
strumental spectral broadening is the same for the tem- 
plate stars and the galaxy spectrum, and we need not 
know the underlying instrumental line profile. 

As possible template stars we chose the follow- 
ing: K3V, MOIII, MOV, M4V, M5II, and M5III. The 
wavelength regions with emission lines were omitted 
from the fit and the optimal t emplate was convolved 
with Gauss-Hermite expansions ([van der Marel fc Franxl 



fl99l lGerhardlH993l) up to hi, minimizing the difference 
between the galaxy spectrum and the template. 

In addition to the stellar continuum, the non-stellar 
continuum is fitted via an additive Legendre polynomial 
of fourth order. This becomes important in the central 
region where the power-law continuum of the AGN dom- 
inates the flux distribution. After the subtraction of the 
stellar and non-stellar continuum we are left with a data 
cube of pure emission line spectra. 

2.3. Extraction of the gas emission lines 

The SINFONI Cen A spectra in the wavelength range 
A1.43-1.87iim and Al.93-2.47/xm (the H- and K-band, re- 
spectively) exhibit a wealth of gas emission lines. High 
excitation lines such as [Sivi] and [CaVlli], ionized gas 
emission lines ([Fen], Bi'7, Hei), and several transitions 
of molecular hydrogen H 2 are detected (Figure [2]) . In 
this paper we focus on the kinematic properties of [Si vi] , 
Br7, [Fen], and H2 and show that they have quite dif- 
ferent kinematics. While the high-excitation lines appear 
to be affected, or created, by Cen A's jet, the H 2 gas ap- 
pears to be solely rotating. We construct a dynamical 
model to explain the kinematics of the strongest line of 
molecular hydrogen 1-0 S(l) H2 at A = 2.122/xm and use 
this to measure the mass of the supermassive black hole 
at the center of NGC 5128. Single Gaussians provide a 
good fit to the emission lines and are used to measure 
the central wavelength, width and intensity of each line 
independently. The fit is performed in IDL 7 using a non- 
linear least squares fit to the line and the errors are the 
l-cr error estimates of the fit parameters. The typical 
uncertainties of the peak flux and the width of the lines 
are 3-5%, whereas the position of the line ca be deter- 
mined to less than 1% accuracy. This translates to a 
typical uncertainty in the velocity of the molecular hy- 
drogen line of ~ ±5 km s _1 and in the velocity dispersion 
of ~ ±10 km s" 1 . 

To allow for the best possible extraction of the gas lines 
the extraction window is centered on the expected wave- 
length. Its width is optimized iteratively to fully cover 
the width of the line and to make sure the extraction win- 
dow covers the same range on the left and right of the 
line peak. An initial estimate of the Gaussian fit parame- 
ters (amplitude, central position, and width) was derived 
from a smoothed spectrum at a central position in the 
velocity field, near the AGN, and applied as a starting 
value throughout the field. This initial estimate is intro- 
duced in order to prevent fitting of spurious lines in the 
regions of bad pixels or strong continuum variation. 

The signal-to-noise ratio of the spectra appears high 
enough for extracting the gas emission lines and drops 
rapidly to zero. Therefore, we do not spatially bin the 
data. We only consider the detection of the lines to be 
secure when their amplitude is a factor of 3 above the 
RMS scatter of the spectrum (A/N>3). In this way we 
get an accurate fit to the lines over the entire field, with 
typical uncertainties of a few percent. 

3. EMISSION-LINE GAS MEASUREMENTS 

From the parameters of the Gaussian fit - peak value, 
mean wavelength and width - we get the flux (or, surface 
brightness), the velocity and the velocity dispersion for 

' See http://www.rsinc.com 
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Fig. 2. — K-band spectrum of the nucleus of Cen A integrated over a square aperture of 0"8 X CK'8. The best-fitting composite spectral 
template (a mix of stellar spectra and a featureless continuum) is overplotted in red. The residual, reflecting the pure emission line 
spectrum, is shown in the lower part, with the main emission lines identified. Note the strong CO absorption that is very well fitted by the 
stellar template 



the considered gas species, [Sivi], Br7, [Fen], and H 2 . 
Since the stellar and non-stellar continua have been sub- 
tracted, the line-flux is directly measured as v / 27rF pGa ]< ; <7. 
For the velocity, we take the recessi on velocity of Cen A's 
stellar body (v sys = 532 ±5 km s~^ lMarconi et alj|200ll) 
as the reference and measure all line shifts with respect 
to this velocity. 

3.1. Gas kinematics 

Figures [3][6] show the maps of total flux, velocity, and 
velocity dispersion for the detected lines of [Sivi], Br7, 
[Fe n] , and H 2 . These line maps illustrate vividly how the 
flux distribution and kinematics change when going from 
high- to low-excitation states. The highest excitation 
line, [Sivi], is dominated by a non-rotational component. 

When comparing the velocity fields of [Sivi], [Fen], 
and H2 , (middle panel in Figures [3][6]) one notices that 
the velocity field of [Si vi] consists of two major compo- 
nents: rotational and translational motion. The velocity 
fields of [Fen] and Br7 are dominated by rotation but 
are still distorted by a non-rotational component that 
is strongest to the lower right (south-west) of the field 
(blue component). These non-rotational motions seen in 
[Si Vi] , [Fe II] , and Br7 are located close to the projected 
direction of the radio jet in Centau r us A (P.A.=51° 
Clarke et al]|1992t iTingav et all 119981 : lHardcastle et all 
20031) . The measured values for the inclination angle of 
the jet vary between 50° < i < 80° (iTingav et al.lll998L 
from VLBI data) and 20° < i < 50° (|Hardcastle et al.1 
120031 from VLA data), but without doubt the North- 
Eastern part is pointed towards us. Surprisingly, the 
direction of motion seen in [Sivi] is red-shifted in the jet 
pointing towards us and blue-shifted in the counter-jet. 
Under the plausible assumption that the geometries of 
the jet and [Sivi], are aligned, this indicates an inflow of 
material towards the nucleus. 

This inflow motion could be associated with backflow 
of gas that was accelerated by Cen A's jet and after pro- 



ducing a bowshock is flowing back at the side of the 
jet cocoon. Although on much large r scales, this phe- 
nomenon is seen in jet simulations (e.g. lKrausdl2005l ) . As 
the gas is flowing back, it is reionized by r adiation from 
the ce ntral source. This was proposed by iTavlor et al.1 
( 1992) as a mechanism to produce the narrow line re- 
gions in Seyfert galaxies. 

While the kinematics of the very high and medium 
ionization lines [Si vi] and [Fe 11] are manifestly influ- 
enced by this jet induced motion, the velocity field of H 2 
shows no distortion due to the jet. The velocity field is 
very smooth and symmetric and indicative of rotation, 
with a striking twist of the major kinematic axis around 
a median P.A.=155°. However, the velocity fields in all 
the other gas tracers give evidence that at least part of 
that gas rotates in a nuclear gas disk. Moreover, the 
velocity dispersion maps of [Sivi], Br7, [Fen], and H2 
support the picture of an inclined nuclear gas disk. They 
all show an elongated structure in their high dispersion 
component, mimicking a disk, at a position angle of 
125 ± 25°, and a declining dispersion profile outwards. 
The ionized gas species ([Sivi], [Fell], and Br7) show 
in addition a colder component (with a ~ 150 kms -1 ) 
that is elongated in the direction of the translational 
motion. Although being as high as a ~ 400 kms -1 , the 
velocity dispersion of H 2 has the lowest central value. It 
is not clear what causes this high velocity dispersion for 
the molecular gas (see also 14. 2[) . Nevertheless, H 2 shows 
the most ordered structure both in the maps of velocity 
dispersion and mean velocity and it seems to be well 
settled in a disk. 

To visualize the difference in the velocity pattern of 
[Si Vi] and H 2 Figure [7] shows the comparison of the H 2 
and [Si Vi] velocity fields (left and middle panel) masked 
with the [Sivi] flux map. The right panel shows the 
direct difference u([Sivi]) — u(H 2 ). Overplotted are the 
VLA contours of the radio jet (unpublished VLA data 
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kindly provided by M. Hardcastle). The non-rotational 
component of the [Si vi] velocity held becomes strongest 
SW of the nucleus, and can be identified with the in- 
nermost knot in the VLA radio jet. We most likely see 
evidence for jet-gas cloud interaction in the ionized gas 
species. 



3.2. Gas morphology 

Looking at the flux maps of the gas species, one no- 
tices the very different morphologies in the high- and low- 
ionization lines. The morphology of [Si vi] is dominated 
by an elongated structure that extends south-west of the 
nucleus, at a position angle of ~ 33°. This structure 
widens and ends in a blob or knot at ~ 1" ~ 17pc from 
the center. The position of this knot is coincident with 



and denoted 



the innermost knot in the radio co unterjet south-west 
of the nucleus (seen bylClarke et al.l (fl9 
SJ1 bv lHardcastle et all (|2003Q ). 

Overall, the morphology of Br7 and [Fe n] (Figures H] 
and [5|) resemble that of [Si vi] very closely, but overall 
the elongation is not as pronounced and the structure 
appears rounder. 

For all gas species, the P.A. of the elongation is ~ 33° 
which is t he same as for the elon gated structure detected 
in Paa by Sch reier et al.l (|1998l ). that is centered on the 
nucleus and extended by ±2". They interpret this as an 
inclined, ^40 pc diameter, thin nuclear disk of ionized 
gas rather than a jet-gas cloud interaction. However, 
our 2D data show that the gas moves along this elongated 
structure that extends around the jet axis. The direction 
of motion hints towards a backflow of jet material onto 
the accretion disk as mentioned above. 

In addition to this elongated structure, our high spa- 
tial resolution integral field data show evidence for a nu- 
clear disk of ionized and molecular gas oriented approxi- 
mately perpendicular to the jet angular momentum vec- 
tor. Looking at the central (X'5 x 0'.'5 of the H2 flux map 
(Fig. 0, one notices a disk-like structure with a major 
axis of ~ 140°. The same structure is visible in the Br7 
and [Fen] flux maps, although a bit rounder. Looking at 
the whole H2 flux distribution, it appears very different 
from the ionized gas species, and the shells at the upper 
left and lower right are a dominant feature. These shells 
are reminiscent of the bowshock structures s een at larger 
scale s in the outer regions of radio jets (e.g. ICarilli et al.1 
119881 for Cy g A) and on smaller s cales in Herbig-Haro 
objects (e.g. iReipurth et all l2002f ). They are located 
at a (projected) position where the elongated structure 
in [Sivi], Br7, and [Fen] disappea rs. This is another 
hint towards the bowshock model of I Taylor et al.l (|1992l ) 
where the shocked gas (here [Si vi] , [Fe 11] , and Br7) flows 
back towards the nuclear source along the shell. It is im- 
portant to note that although the shells dominate the 
appearance of the H2 flux map, they do not leave any 
kinematic signature in the H2 velocity field. 



4. GAS DYNAMICAL MODELING 

While the (highly) ionized gas shows (strong) kine- 
matic influence by the jet, the molecular gas (H2) seems 
to "feel" only gravity. This is the reason why we focus 
on H2 when we construct a dynamical model to measure 
the mass of the central supermassive black hole. 



4.1. Method 

The dynamical model follows the approach of HN+06 
and uses two-dimensional SINFONI gas kinematic maps 
as constraints. In brief: to explain the H2 gas mo- 
tions seen in the center of Cen A we construct a kine- 
matic model where we assume the gas moves in a thin 
disk solely under the gravitational influence of the sur- 
rounding stars and the expected central black hole. 
Then, the gravitational potential is given as $(r) = 
$*(r) + ^bhM- The stellar potential is taken from 
HN+06, where NaCo, NICMOS, and 2MASS K-band im- 
ages of NGC 5128 are used to construct a Multi-Gaussian 
Expansion (MGE) para meterization to the surface 
brigh t ness of this galax y (fEmselle m. Monnet. fe; Baconl 
Il994t ICappellarll l2002al ). The assumptions of spher- 
ical symmetry and constant stellar mass-to- light ra - 
tio (M/L K = (0.72 ± 0.04) M /L dSike et al.l 1200^ : 
HN+06) lead to the three-dimensional mass model, that 
gives the stellar velocity contribution to the dynamical 
model. 

Our dynamical model is based on the widely 
used approach to model the emission line profile of 
gas moving m a thin disk dMacchetto et all 119971; 



van der Marel fe van den Boschl [1998E Bertola et al.l 
19981: iBarth et al.ll2001l ). In HN+06 we considered three 



conceptual modeling approaches: 1) a cold disk model 
that fully neglects the velocity dispersion, 2) a hot disk 
model that accounts for the high velocity dispersion of 
the gas, and 3) a spherical Jeans model that accounts for 
the high velocity dispersion but neglects the indicated 
disk geometry. We do not repeat the spherical Jeans 
model here, as the H2 velocity maps clearly require a 
disk model. Since the observed velocity dispersion even 
of the H2 gas at the nucleus of Cen A exceeds the mean 
rotation by more than a factor of two, we must account 
for the velocity dispersion in the dynamical model. We 
assume the gas disk to be geometrically flat but with an 
isotropic pressure and construct an axisymmetric Jeans 
model in hydrostatic equilibrium. In this case, the mean 
rotation velocit y (azimuthal velocity) v$ i s given by the 
Jeans equation (|Binnev fe Trem aine 1987t, Eq.4-64a) 



R d{p g (j 2 R ) 



v<;> = 



Pa 



OR 



R 



dR 1 



(2) 



where R is the projected radius and an is the radial 
velocity dispersion of the gas. We assume that the H2 
surface brightness of the gas disk E g reflects the tracer 
gas density p g (see Section 14.31 for a discussion) . Note 
that the model is not self-consistent, i.e. the contribu- 
tion of the gas mass to the overall potential is neglected. 
We estim ate the mas s insid e the SINFONI field-of-view, 
following Isra el et "ail (jl990) . They identify a central H2 
gas disk with an outer radius of ~ 6'.' 5, a thickness of 
~ 3'.'3, and an inner cavity of ~ l'/7 (assuming a distance 
to Cen A of D=3.5 Mpc), and with a density distribution 
of p(r) ~ r~ 2 . The total disk mass is 2 x 10 7 M o . To 
get a first order estimate of the mass inside the central 
3" x 3" we assume a constant density throughout the disk 
and neglect the inner cavity, which is not confirmed by 
our data. Inside the field of view of our SINFONI obser- 
vations we therefore expect Mjj 2 ~ 1 x 10 6 M Q , which 
is a factor of 60 to 200 smaller tha n the black hole mas s 
measurements for Cen A (HN+06: iMarconi et al.ll2001h . 
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Fig. 3. — Surface brightness, velocity, and velocity dispersion maps of the [Sivi] line emission (left, middle, right panel, respectively). 
The surface brightness map is displayed in logarithmic scaling. The velocities and velocity dispersions are only plotted for the regions 
that are above the lower flux limit displayed in the left panel. The velocity scale is given relative to the systemic velocity of the galaxy 
(532 km s" 1 ). The velocity dispersion is corrected for the instrumental dispersion (<Tj nstr ~ 65 km s — 1 ). North is up and East is to the left 
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Fig. 4. — Flux, velocity, and velocity dispersion maps of Br7. See caption of Fig. [3] for more details 
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and therefore negligible. 

To match the observations, the resulting velocity field 
can then be projected onto the plane of the sky, given 
the inclination angle i and the position angle (P. A.) of 
the projected major axis of the gas disk. For compar- 
ison with the data, the projected velocity field must 
also be broadened with the velocity dispersion and then 
weighted by the gas surface brightness, using the param- 
eterizations given subsequently in Sections 14.21 and 14.31 
respectively. Finally, we simulate observations of the disk 
through the SINFONI instrument, i.e. we convolve with 
the SINFONI PSF and sample over the pixel size, to 
achieve the best possible match to the data. The model 
thus has three free parameters, the black hole mass Mbh, 
the inclination angle of the disk i, and its projected P.A. 
(. The galaxy center position on the detector (defined 
as the continuum peak) and the stellar mass-to-light ra- 
tio are fixed beforehand, with M/Lk = 0.72 M Q /L0 as 
described above. To measure the black hole mass in 
NGC 5128, we run models with different values for the 
free parameter and look for the best possible match to 
the data, minimizing the 2. 

For the modeling we use the IDL software 8 of HN+06, 
which accounts for the SINFONI PSF (as determined 
in Section l2.ip . instrumental broadening, and the finite 
SINFONI pixel size, to generate a two-dimensional model 
spectrum with the same pixel scale as the observations. 

8 using Craig B. Markwardt's MPFIT package 



The extraction of the mean velocity and velocity dis- 
persion from this synthetic data cube is carried out in 
exactly the same manner as for the observed data, by 
fitting single Gaussians to the individual spectra at each 
pixel. 



4.2. Intrinsic gas velocity dispersion 

The observed velocity dispersion of the molecular hy- 
drogen H2 gas at the center of NGC 5128 peaks at 
400 km s _1 and exceeds the measured mean rotational 
velocity by more than a factor of two (for compari- 
son: the ionized gas species [Fen] used in the study of 
HN+06, has a peak velocity dispersion of ~ 600 km s _1 ). 
The physical origin of this high velocity dispersion is 
not clear. Partially it might be expla i ned by spa - 
tially unresolved rotation (|Marconi et alj I2001L l2006j) . 
or it might be due to local turbulent gas motions, as 
suggested by several authors for other galaxies (e.g . 
[ Barth et aljfeOQlt Ivan der Marel & van den Boschlll998l; 
IVerdoes Kleiin et al.ll2002f ). which appears problematic 
especially for the H2. 

Regardless of its physical cause, we consider the high 
velocity dispersion to contribute to the pressure support 
of the gas disk. Following the approach of HN+06 we in- 
clude the velocity dispersion to the gas dynamical model 
via an isotropic pressure term in the Jeans equation. The 
rotational velocity therefore becomes sub-Keplerian. We 
find that the intrinsic velocity dispersion is well described 
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Fig. 5. — Flux, velocity, and velocity dispersion maps of [Fell]. See caption of Fig.[3]for more details. We masked out an artifact around 
y~ — CK'7 that is caused by illumination effects on the SINFONI detector 
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by a double exponential profile of the form 



OR = cr e 



—r/ri 



(3) 



and we fit the observed H2, velocity dispersion profile 
for the best set of parameters (00 = 140 km s -1 , o\ = 
65 km s _1 , ro — 0'.'25, and r\ = 4'.'0) to get the intrinsic 
dispersion profile. 

4.3. Emission line surface brightness 

For various aspects of the modeling, we need to know 
the intrinsic spatial emission line profile (assumed to be 
axisymmetric) . But a direct deconvolution of the ob- 
served surface brightness is very difficult, since for adap- 
tive optics observations the exact shape of the PSF is 
unknown and we would therefore introduce artifacts to 
the gas distribution that bias the velocity distribution 
inside the inner 0'.'5. 

Furthermore, the gas morphology of H2 is quite com- 
plex and it is not a priori clear how this complex morphol- 
ogy leads to such a smooth velocity field. On the other 
hand, the surface brightness does not necessarily resem- 
ble the real physical gas structure. The bulk mass in H2 is 
most probably at very low temp eratures (T~ 10 — 15 K, 
as derived by llsrael et al.l 1 19901 from CO observations) 
and therefore not excited. 

Figure [5] shows slices through the data cube around the 
H 2 line. The width of the velocity slices is 33 km s _1 , 
corresponding to one pixel in wavelength direction (i.e. 
half the spectral resolution cxinstr)- Material located near 



the nucleus is present in several consecutive panels, indi- 
cating that the material has a high velocity dispersion. 
In addition there is material that appears in shell-like 
structures only in one or two bins, hinting towards a 
lower velocity dispersion. The surface brightness of the 
H2 gas seems to be highest along the jet direction for 
0'.'5 < r < l'/3 (see Figure EJ). We actually see shell- 
like structures that might be due to shocked H2 gas; this 
hypothesis is supported by the fact that the [Sivi] and 
Br7 gas distributions fit quite nicely into the shells seen 
in H 2 . 

Looking at the flux distribution of H2 (Figure [6]) , in 
the radial range r < ff.'i a elongated structure which 
might be reminiscent of a disk is visible with a major 
axis of P.A.=136°. Ellipse fits to this disk structure give 
a minor-to-major axis ratio of q ~ 0.67 which translates 
to an inclination angle of ~ 48° given a circular thin 
disk configuration. The hypothesis of the central disk 
structure is also supported by the shape of the velocity 
dispersion field. Looking at the right panel of Figure [5] 
we see a disk-like structure with a major-axis position 
angle of ~ 140° inside ~ 0'.'4. 

Drawing on this qualitative description of the H2 line 
distribution, and the conjecture that the gas moves in a 
disk, we model the emission line surface brightness in two 
parts: an exponential disk that dominates the inner 0'.'5 
and a smoothed version of the actually observed H2 flux 
distribution resembling the detailed gas morphology in 
the region outside 0'.'5, where the PSF convolution is less 
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Fig. 7. — Comparison between the [Sivi] (left) and H2 (middle) velocity fields with their difference shown in the right panel. The 
velocity fields are masked with the [Si vi] flux map which nearly coincides with the radio jet (red contour overlay [unpublished VLA data 
kindly provided by M. Hardcastle]). The differenced velocity field (right hand panel), meant to 'correct' the [Sivi] kinematics for their 
rotational component, reveals a distinct [Sivi] velocity gradient along the jet. The blue component SW of the AGN (0'.'7,— 0'.'6) marks the 
strongest translational velocity component, matching the knot in the radio jet 



critical. For the black hole mass modeling it is not im- 
portant to add the second (larger) component, since the 
inner < 075 are dominant. Nevertheless we add it to get 
a better estimate of the influence of the non-symmetric 
gas distribution on the appearance of the gas velocity 
field, since part of the small-scale structure in the ve- 
locity fields can be reproduced b y a patchy gas su rface 
brightness folded into the model ( Barth et al.l 200ll ). 

We have tested the influence of the parameterization 
of the exponential disk component on the resulting best- 
fit black hole mass, and find that a profile that puts 
^12% more (less) flux inside the central 0'.'5 (but still 
fits the observed flux distribution well) results in a black 
hole mass that is less than 3% lower (higher). This re- 
sult is in line with the extensive tests on the influence 
of the surface br i ghtnes s parameterization carried out 
bv iMarconi et all |2006). In any case, this uncertainty 
is small compared to other uncertainties that enter the 
black hole mass measurement, e.g. the inclination angle 
of the gas disk. 

4.4. Tilted-ring model 

The kinematics predicted by a flat, or co-planar, thin- 
disk model must be an oversimplification, as becomes 
apparent when comparing to the observed data. The 
twists in the velocity field cannot be reproduced, and the 
velocity gradient is strictly declining from the peak at ~ 
C//2 outwards without being able to resemble the second 
and third peak at 0'.'45 and l'.'f (~ f9pc). We therefore 
mode l the kinematics via a tilted-ring model (Begeman 
119871) as it was done b efor e on scales of r ~ 0.'5 to 3' 
bv lQuillen et alj (| 19921 ) and iNicholson et all (| 19921) who 
modeled the CO (2-1) and Ha velocity fields in Cen A, 
respectively. The difference to a coplanar model is that 
the inclination angle and position angle of the gas disk 
are a function of radius. The orbits of the gas at each 
radius remain circular, but neighboring orbits are not 
necessarily in the same plane. The gas-disk geometry 
changes from co-planar to warped. 

We work in polar coordinates and use discrete radial 
steps where the model is to be calculated. The model is 
linearly interpolated between the discrete points on the 
model grid. The gas disk is made up of concentric rings. 
Each ring is represented by three parameters: its radius, 
R, inclination angle i, and azimuthal angle £ (relative 
to the projected major axis). When projected along the 
line-of-sight, the rings become ellipses. The flattening, q, 



of the major and minor axes is related to the inclination 
angle i via cos(i) = q. The flattening defines an ellipse 
on the sky of ellipticity e = 1 — q. If the gas is assumed 
to move on circular orbits along the rings, the projected 
velocity can be described by the simple cosine form 

v(R, C) = v sys + v c (R) sm(i(R)) cos(((R)), (4) 

where v c (R) denotes the circular velocity on a given ring 
of radius R, while v sys gives the systemic velocity of the 
entire galaxy. This method of expanding the full gas ve- 
lo city field in a set of tilted rings goes back to the work 
of iBegemanl (Il987[). H ere, we are using the method of 
IKrainovic et alJ ([2006) , called kinemetry, to determine 
the set of ellipses along which the velocity field is best de- 
scribed by equation ^ . This method is a generalisation 
of surface photometry to all moments of the line-of-sight 
velocity distribution, ft performs harmonic expansion of 
2D maps of observed kinematic moments (velocity, veloc- 
ity dispersion, and higher Gauss-Hermite terms) along 
the best fitting ellipses (either fixed or free to change 
along the radii). The parameters of the best-fit kineme- 
try model, using 27 tilted rings, are plotted in Figure [9] 
and listed in Table [U while Figure [10] shows the fitted 
circular velocity map, using these parameters. Note that 
the parameters of the tilted-ring model vary smoothly 
with radius, although we do not restrict the kinemetry 
routine to smooth functio ns i(R) and PA(R). W e re- 
fer the interested reader to IKrainovic et a l. (|2006f ) for a 
detailed description of the kinemetry method. 

5. RESULTS 

We now focus on modeling the H2 gas kinematics in 
the central 1" (~ 17pc) of NGC 5128 with a warped 
disk model where the gas orbits in the joint potential of 
the stars and the central black hole. The central black 
hole mass as well as the overall inclination angle of the 
gas disk are varied to find the best-fit to the observed 
velocity and velocity dispersion maps. 

5.1. Structure of the H2 disk 

The rotational velocity field seen in molecular hydro- 
gen shows beautiful symmetry about the center despite 
the twist that alters the projected orientation of the kine- 
matic major axis by ±14° about a mean value of ~ 155° 
(see TablcQ]). Especially for the central values, the orien- 
tation of the H2 gas disk is consistent with an orthogonal 
disk-jet picture, as th e radio jet in Cen A is located at a 
position angle of 51° (jTingav et al.lll998f) . 
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Fig. 8. — Velocity channel maps for the molecular gas (H2 ) at A2.121 /an observed in the K-band. The middle panel shows the flux in 
the velocity bin that corresponds to the recession velocity of the overall galaxy (532 km s — 1 ). The other bins slice the cube at blue- and 
redshifted velocities to this zero position. The width of the velocity slices is 33 km s — 1 



In addition to the twist in position angle there is a 
superimposed variation of inclination angle that causes 
the maximum of the projected rotational velocity to vary 
by ± 30 km s _1 (around a mean value of 90 km s _1 , cf. 
Table [T]) . This variation is indeed significant compared 
to the uncertainty of the velocity which is ~ ±5 km s _1 
(see Section POl 

The structure of the gas disk as fitted by the tilted- 



ring model is presented in Figure 1111 For this fit the 
velocity field was symmetrized using the assumption 
v(x,y) = —v{—x,—y). The left panel (a) shows the ac- 
tual kinemetry fit to the symmetrized velocity field. In 
panel (b) the disk is rotated by 45° to an almost edge- 
on view, to make the warp clearly visible. The position 
angle twist is indicated by the red and blue lines for the 
receding and approaching side, respectively. 
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Fig. 9. — Disk parameters of the kinemetry fit to the velocity 
field as a function of projected radius. The inclination angle i is 
derived from the flattening q (i = arccos(g)) 

Our SINFONI integral field data allow us to constrain 
the inclination angle of the gas disk very well. Its mean 
value is 45° ± 7°, which is in very good agreement with 
the best-fitting inclination angle measured by HN+06 
from their NaCo long-slit data (45°). The inclination 
angle is also in agreem ent with the value derived by 
lHardcastle et al.1 (|2003ft from VLA data (20° < i < 50°), 
it is however so mewhat smalle r than the value from VLBI 
data derived by iTingav et~aT1 (j!998f ) (50° < i < 80°). 

The central H2 gas kinematics are well described by 
the tilted-ring model. As already stated above, a copla- 
nar disk model would not be able to reproduce the twist 
in the rotational velocity field and would therefore rep- 
resent an oversimplification of the kinematic structure. 
The coplanar disk model can be ruled out at almost 3<7 
confidence level. 

The warp of the larger scale disk (> 10") was probably 
caused by the merger event, that occurred a few times 
10 s years ago (jQuillen et al.lll993t |Pe"ng et al.ll2002fl . It 
is not a priori clear whether the warp in the innermost 
arcsecond is connected to this larger scale warp that cre- 
ates the prominent a ppearance of the dust disk in Cen A 
(Quill en~et al.|[2006l ). For the nuclear gas disk on scales 
5, 1"2 (< 20pc ) self-induced warping of the accretion disk 
(|Pringlej Il996f ) might be the driving mechanism. Our 
model is not intended to explain the origin of the warp, 



TABLE 1 

Parameters of the tilted ring model 



an 


p.a. [°; 


q 


i[°] V 


[km s" 1 ] 


0.05 


144.0 


0.71 


45.0 


55. 


0.10 


148.5 


0.77 


37.6 


106. 


0.16 


152.1 


0.79 


38.8 


112. 


0.22 


154.8 


0.70 


42.2 


107. 


0.27 


160.5 


0.66 


44.5 


110. 


0.33 


165.4 


0.56 


42.3 


110. 


0.39 


167.8 


0.60 


44.1 


113. 


0.45 


168.5 


0.64 


45.7 


115. 


0.51 


168.9 


0.67 


44.0 


109. 


0.57 


169.0 


0.78 


39.5 


103. 


0.63 


167.2 


0.81 


36.1 


96. 


0.69 


165.2 


0.82 


35.1 


89. 


0.76 


163.5 


0.82 


34.2 


82. 


0.82 


161.7 


0.82 


35.2 


79. 


0.89 


159.8 


0.80 


36.9 


76. 


0.96 


157.0 


0.80 


36.8 


76. 


1.03 


153.8 


0.80 


36.2 


75. 


1.10 


151.1 


0.80 


36.8 


75. 


1.18 


149.3 


0.76 


39.7 


76. 


1.26 


147.8 


0.74 


42.5 


77. 


1.34 


146.8 


0.69 


46.7 


80. 


1.42 


145.1 


0.65 


48.9 


83. 


1.51 


144.0 


0.61 


51.6 


85. 


1.60 


142.8 


0.58 


53.6 


87. 


1.69 


141.6 


0.56 


55.6 


88. 


1.79 


137.9 


0.64 


54.6 


90. 


1.90 


134.9 


0.58 


58.9 


89. 



but is rather meant as a geometric model to resemble the 
observed velocity field. 

5.2. Importance of the inclination angle 

The best-fit black hole mass in gas dynamical models 
depends strongly on the inclination angle assumed for 
the gas disk model (Mbh ^ sin(i)). Any uncertainty in 
the disk inclination angle will directly propagate to the 
uncertainty in the black hole mass. 

Our tilted-ring model fitted by kinemetry fixes the 
overall shape of the gas disk in the model. However, 
we allow for overall variations in the inclination angle, 
not to overconstrain the model. That means the posi- 
tion angles for the tilted-rings remain fixed, as well as 
their relative orientation. While we fix the overall 'twist' 
of the disk a priori, we allow the full dynamical model 
to re-fit the global inclination. In that way we make use 
of the geometrical information that is contained in the 
velocity field (using kinemetry), but at the same time 
give our dynamical model the freedom to find the best 
inclination angle given the general assumptions of the 
model. 

5.3. Best-fit model and black hole mass 

We calculate a grid of possible models for varying disk 
inclination and central black hole mass to get the set 
of values that best match the observed data. The best- 
fitting black hole mass in our tilted-ring model to the 
H2 kinematics is Mbh = x 10 7 M q for a median 

inclination of ~ 34° ± 4° (error bars are given at the 
3(7 level). The best model has a 2 m ; n of 8.2. This rep- 
resents the minimum in the A2 distribution, shown in 
Figure 1121 The points represent models. The contours 
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Fig. 10. — The symmetrized H2 velocity field (left panel) compared to the best fit kinemetry model (right panel). See Section 14.41 for a 
description of the kinemetry model. The ellipse shows the flattening of the largest circle 
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Fig. 11. — Tilted-ring model to describe the nuclear H2, gas disk seen at (a) the original orientation as fitted by kinemetry and (b) rotated 
by 45° to make the warp more visible. The red and blue lines are the line-of-nodes for the receding and approaching side, respectively 



were determined by a two-dimensional smoothing spline 
interpolated from these models and represent A2 values 
of 1.0, 4.0, and 9.0. This corresponds to 68.3%' 95.4%, 
and 99.7% confidence levels for 1 degree of freedom, or 
1(7, 2(7, and 3c confidence levels, respectively. The as- 
sociated best-fit model velocity maps are shown in com- 
parison to the data in Figure [T31 If we would keep the 
inclination angle fixed to the mean value of 45° given by 
kinemetry the model is not able to reproduce the overall 
velocity field. The main deviation between model and 
data appears for radii larger than ~ 0'.'7, i.e. outside the 
radius of influence of the black hole, where the model 
results in a rotational velocity that is significantly higher 
than the observed one. The 2 m i n for this "best-fit" model 
is 16.8 with a best-fit M B h ~ 3.0 x 1O 7 M . This model 
is thus ruled out to over 3er. The entire velocity field 
is better reproduced by a disk at lower inclination (i.e. 
more face-on) plus a higher black hole mass that makes 
up for the decrease in velocity inside the radius of influ- 
ence of the black hole. 

Figure [TJ] shows a comparison of the model and the 
data for the case of no central point mass. Here, the 
gravitational potential is made up only by the stars. The 



mass-to-light ratio is O.72M/L 0iK (HN+06). It is obvi- 
ous that this is no good fit in the central 0'.'5 x 0'.'5. The 
modeled rotation only catches up with the data outside 
~ l'.'O, where the stars clearly dominate the gravitational 
potential. The case for no black hole is excluded to very 
high significance (over 8a). 

There are various factors that influence the black hole 
mass estimate in our dynamical model and we have done 
a substantial number of tests to scrutinize their impact 
on the best-fit result. As mentioned in Section 15.11 the 
assumed geometry of the disk (warped vs. flat) has a 
small influence on the black hole mass. The same holds 
true for the parameterization of the surface brightness. 
We modeled the kinematics for three different parame- 
terizations of the disk's surface brightness profile, with 
all three being a reasonable fit to the data. We found 
that the black hole mass does change by less than 3% 
depending on the assumed surface brightness profile of 
the inner gas disk. This result is in agre ement with the 
detailed analysis of M arconi et a l. (2006). 

Obviously, also the contribution of the stellar poten- 
tial to the total gravitational potential influences the 
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resulting best-fit black hole mass. We used the two 
ext reme value s 0.72M /L w .k and 0.53M/Lq,k derived 
by ISilge et al] (|2005h through stellar dynamical mod- 
els at inclination angles 90° and 45°, respectively (the 
value for 20° is 0.68M/L Qj k and lies in between these 
two). Using 0.53M/L Qi k, the best-fitting black hole 
mass increases by ~ 8% compared to the best-fit value 
of M BH = 4.5 x 10 7 M Q (for 0.72 M/L Q , K ). The uncer- 
tainty introduced by the inclination angle is much larger 
than this, which is reflected by the fairly large uncer- 
tainty limits on our black hole mass measurement. 

5.4. Asymmetries in the H2 velocity field 

The overall shape of the H2 velocity field appears to 
be point-symmetric about the position of the AGN (the 
peak in the K-band continuum which is coincident with 
the peak in the H2 surface brightness). This symmetry 
is proven and illustrated by the smooth kincmetry fit. 

However, taking a closer look at the peak velocities in 
the field, it is striking, that the peak velocity in the blue 
dip at (x,y) = (— 0'.'2, — 0'.'5) exceeds the one in the cor- 
responding red peak ((x, y) — (0'.'2, ft' 5)) by ~ 30 km s" 1 
(corresponding to ~ 30% of the peak velocity) . 

The reason for this asymmetry is not merely an over- 
all velocity shift, since that would affect the whole field. 
This can be ruled out by the fact that the absolute val- 
ues of the velocities are in excellent agreement for radii 
larger than ~ ft'8. 

It is instructive to overplot the iso-flux contours to the 
velocity field (see Figure H~5]). since one clearly sees that 
the asymmetry in the velocity field is also present in the 
flux distribution. This leads us to speculate that the 
gas density is not the same throughout the disk, as a 
higher density leads to a higher flux. The density might 
be increased by gas that is streaming towards the central 
disk. This might be of similar origin as the non-rotational 
motions observed in [Fen], Br7, and [Sivi] but is not 
as obvious in molecular hydrogen. The asymmetry is 
also present in the Br7 and [Sivi] velocity maps, being 
strongest in [Si vi] . What we observe here, might be the 
fueling process of the nuclear disk. 

Our gas dynamical model is by construction symmet- 
ric, and therefore not able to reproduce asymmetries 
in the velocity structure. For determining the best-fit 
black hole mass, we have minimized the difference of 
the symmetric model to the asymmetric data, looking 
for the best compromise. After having gained the best- 
fit model (to the original data) we are free to compare 
this to a symmetrized version of the H2 velocity field, 
v sym (x, y) = (v(x, y) - v(-x, -y))/2. The agreement, as 
shown in Figure \16\ is astonishing. 

6. DISCUSSION 

The high resolution 2D data have now reached a qual- 
ity that the existence of a central black hole (or compact 
mass) can be demonstrated, and its mass be estimated, 
even without a sophisticated model, e.g. by just using 
the concept of the radius of influence of the black hole, 
tbh- This is the region where the mass of the enclosed 
stars equals the mass of the supermassive black hole: 



where Mbh is the mass of the black hole, G is the grav- 
itational constant and ex* is the velocity dispersion of 



the stellar spheroid. In the velocity field, this is the 
point of minimum rotation, where the black hole stops 
to dominate the gravitational potential, before the stars 
(with their rising rotation curve) take over. For our 
H2 velocity this is at ~ 0'.'8 as seen in Figure [16] (left 
panel). For Cen A we have the f ollowing numbers: 
cr* = 138 km s" 1 (jSilge et al.ll2005l ). at D=3.5 Mpc, 
i'bh = ft'8 ± O'.'l = 13.2 pc±1.7 pc, and therefore we 
get M BH = (6.0 ± 0.7) x 10 7 M . 

The observed radius of minimum rotation is inde- 
pendent of the inclination and so this simple concept 
provides a nice check on the black hole mass derived 
via dynamical modeling. Given the excellent spatial 
resolution of our AO-assisted data (FWHM corc = O'.'l 2 
and FWHMj,„i„ = 0'.'30) the observed radius of influence 
of the black hole is well beyond the radius where PSF 
effects start significantly affecting the derived velocity 
curve. 

The best-fit black hole mass derived through modeling 
of the H 2 kinematics, M B h = (4.5±^o) x 10 7 M Q at 
i = 34° ±4°, is in good agreement with the mass derived 
by HN+06 ((6-l±o;|) x 1O 7 M at i = 45°) using high 
spatial resolution kinematics of [Fe 11] derived from 
AO-assisted NaCo long slit data (FWHM= ft' 11). The 
dynamical model they used is in principle identical 
to the one described above, except the fact that we 
cover the velocity field in two dimensions and they 
modeled only four slit positions. Moreover, they did not 
include the disk inclination angle as a free parameter, 
which is the main reason for their smaller error bars. 
Concerning the disk geometry, HN+06 excluded disk 
inclination angles below 45° d ue to the jet inclination 
derived bv lTingav et~aTl (| 1998ft . who gi ve 50° < i < 80°. 
Howev er, following the analysis of lHardcastle et al.l 
(2003), who use jet-counterjet ratios and apparent 
motions in the jet, the jet inclination is most likely in 
the range 20° < i < 50°. HN+06 state that if they were 
to allow an inclination angle of i — 25° with respect 
to the line-of-sight, their best-fit black hole mass is 
~ 1.2 x 10 s Mq. This is significantly larger than the 
value we measure. But one has to take into account, 
that the modeled gas species are not the same for the 
two studies. While the molecular gas has a well-ordered 
rotation field, the 2D velocity field of [Fen] (Fig. 
clearly exhibits two components: inflow and rotation. 
This superposition could not be seen in the long-slit 
data and was not accounted for in the model of HN+06. 
They modeled the total [Fen] kinematics under the 
assumption of gas rotating in a flat, thin disk. 

6.1. Comparison to previous gas dynamical models 

The agreemen t of our modeling resu lts with the re- 
cent analysis of iKrajnovic et alJ (|2007f ) is comfortable: 
they modeled integral-field Pa/3 kinematics and derived 
M BH = 8.251^25 x 10 7 M Q at i = 25° (where the error 
bars are also at the 3<t confidence level). The resolution 
of their seeing limited data is 0'.'5 — 0'.'6, i.e. a factor of 
~ 4 — 5 larger than ours, but good enough to resolve the 
radius of influence of the black hole in Cen A. Again, 
the somewhat larger black hole mass value migh t be at- 
tributed to the fact that (jKrainovic et al.ll2007ft model 
the kinematics of the ionized gas species Pa/3, which 
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Fig. 12. — Constraining the mass of the central black hole: The Figure indicates the grid of models (in black hole mass, Mgn, and 
disk inclination) that was calculated, and the contours show A2 in the vicinity of the best fit dynamical models for matching the H2, 
kinematics. The minimum 2 model is at a Mbh ~ 4.5 X 10 7 Mq and a median disk inclination of 34°. The contours indicate the la, 2a, 
and 3a confidence levels, respectively (see text for details) 
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Fig. 13. — Velocity field of the best-fitting dynamical model (center, top), with a black hole mass of 4.5 X 1O 7 M0 and a median disk 
inclination of 34°, in comparison to the data (top left). The velocity residual (data-model) is shown in the right panel. The velocity curves 
in the bottom panels are extracted along the linc-of- nodes (ovcrplotted to the velocity maps), and represent the peak velocity curves. 
The diamonds correspond to the model velocity curve while the crosses correspond to the data. A mismatch in the data and model most 
prominent beyond r< — 1'.'7 is likely due to the fact that the inclination angle of the modeled gas disk is not well represented in the 
outermost parts. This has no impact on the derived black hole mass 
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Fig. 14. — Kinematic evidence for a central black hole: The top panels compare the observed symmetrized H2 velocity field (left) to a 
velocity field model for the case of zero black hole mass with only a stellar potential, derived by HN+06 (right). The bottom panels show 
the velocity curves extracted along the line-of-nodes (indicated in the top panels). Obviously, this model is not a good fit to the data in 
the central 0'.'8 X CK'8. Only beyond l'.'O are the model velocity curve (diamonds) in reasonable agreement with the data (crosses) 
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Fig. 15. — The H2 surface brightness contours are overplotted 
to the H2 velocity map. The velocity field is asymmetric in the 
sense that the blue dip is of larger amplitude (in the center-of- 
mass system) than the red peak. This asymmetry is also seen in 
the flux contours, but such asymmetries are not captured by our 
model 

might also be affected by non-gravitational motions. 

This decrease in black hole mass compared 
to recent dy namical mea su remen ts of ion- 
ized gas by IMarconi et"aH (|2001l) (M B h = 

2.0±5;2 x 1O 8 M ), lHaring-Neumaver et all (poM ) 
(M BH = 6.1±g;| x 10 7 M Q ), IMarconi et all (f2006h 
(M BH = 6.5i°;^ x 10 7 M Q ), and iKrainovic et alj (|2007f ) 
(M BH = 8.25^4 25 x 10 7 M Q ) bring s Centaurus A in 
agreement with the Mbh-ct relation (|Ferrarese fc Merrittj 
120001 iGebhardt etldl l2000f) (see Fig. [Tp. All previous 
gas dynamical studies were based on ionized gas species, 
which are most probably all influenced by Cen A's jet, 
and are no good tracers of the gravitational potential in 
the innermost arcsecond. One of the main advantages of 
our high resolution two-dimensional data is to be able to 
reveal the jet's influence and to choose the appropriate 
gas species for the dynamical model. 

6.2. Comparison to stellar dynamical models 

Using stellar kinematics (v, er, /13, and /14) derived from 
Gemini GNIRS data along two sli t-positions with see ing- 
limited resolution of 0'.'4 to 0'.'65 ISilge et alj pOOBT ) de- 
rive a best-fit black hole mass of 2.4+q'J x 10 8 M Q for an 
edge-on model, 1.8to!4 x 10 8 M Q for a model with i=45°, 
and 1.5^Q 2 x 10 8 M© for a model with an inclination of 
20° (note that theses error bars are at the la confidence 
level) . 

Using integral- field stellar kinematics (v, a, /13, and 
/14) derived from the same SINFONI data set used for 
the underlying gas dynamical study, Cappellari et al. (in 
prep.) get a best-fit black hole mass of 5 ± 3 x 10 7 M Q 
(3cr error bars). This value is mo re than a factor o f ~ 
3 lower than the value derived bv ISilge et al.l (|2005l ). it 
is, however, in full agreement with the gas dynamical 
measurement presented in the current paper. 

This is an important result for the general compar- 
ison of stellar and gas kinematical modeling of black 
hole masses. Up to now only a couple of these com- 
parative cases have been investi gated and in half of the 
cases with negative results (e.g. ICappellari et alj|2002bt 
IShapiro et al.ll2006[) . 

This is the first case where the modeling techniques can 



be tested on stellar and gas kinematical data extracted 
from the same complete data set. The fact that we deal 
with integral-field kinematics allows the constraints on 
the black hole mass to be significantly tightened com- 
pared to long-slit observations. The excellent agreement 
between gas and stellar dynamical modeling results thus 
provides stringent evidence for the reliability of both 
techniques. 

7. CONCLUSIONS 

This work presents observations of the nearby active 
elliptical galaxy NGC 5128 (Cen A) with the adaptive- 
optics assisted integral field spectrograph SINFONI at 
the VLT. Our K-band data used to measure the black 
hole mass in Cen A from emission line kinematics have a 
spatial resolution of 0'/12 and an estimated Strehl ratio 
of - 17%. The field of view is 3" x 3". 
In our H- and K-band data we detect the following emis- 
sion lines: [Fell], [Sivi], He I, Br7, [Cavin], and several 
vibrational transitions of molecular hydrogen, H2 (the 
strongest is 1-0 S(l) at 2.12/im). The main results of our 
analysis are: 

• Analysing the velocity fields of [Si Vl] , Bi'7, [Fe 11] , 
and the strongest H2 line, we find that the surface 
brightness and also the motion of the gas species 
is increasingly influenced (or produced) by the jet 
when going from low to high excitation lines. The 
velocity fields of [Sivi], Br7, and [Fen] clearly ex- 
hibit two components: 1) rotational motion at a 
major angle of ~ 150°, consistent with an orthog- 
onal disk-jet picture, and 2) non-rotational motion 
along the direction of the jet, that is consistent with 
a back-flow of gas along the side of the jet's co- 
coon. This non-rotational component is strongest 
for [Sivi]. 

• The surface brightness of H2 shows an indication 
of a central gas disk at a position angle of P.A.~ 
136° and a minor-to-major axis ratio of 0.67. The 
bright spots NE and SW of the nuclear disk remind 
of the lobes seen on much larger scales in jet-gas 
interactions. 

• The overall velocity field of H2 shows beauti- 
ful point- symmetry about the unresolved nucleus. 
However, the H2 peak velocities are asymmetric 
inside the central Of! 7. This asymmetry is also seen 
in the surface brightness map of H2 , and might be 
explained by denser gas that leads to an increase 
in flux. The density could be increased by gas that 
is streaming towards the central disk. This might 
be similar to the non-rotational motions observed 
in [Fen], Br7, and [Sivi] but is not as obvious in 
molecular hydrogen. It is possible that we see the 
fueling process of the nuclear disk. 

Another influence of the jet on the velocity field 
might be the twists in the zero velocity curve (at 
the positions where the lobes appear in the surface 
brightness map of H2). The rotational major axis 
(line-of-nodes) runs from SE to NW, with the SE 
side blue- and the NW side redshifted with respect 
to the nucleus. It follows the shape of an "S" . The 
velocity field is reminiscent of a warped disk. Judg- 
ing from the smooth velocity field of the H2 gas, 
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Fig. 16. — The symmetrized velocity map (top left) in comparison to the model (Mbh = 4.5 X 1O 7 M0 and < % >= 34° - top middle 
panel). The residual (data-model) is shown in the right panel. The velocity curves in the bottom panels are extracted along the line-of- nodes 
(overplotted to the velocity maps), and therefore represent the peak velocity curves. The diamonds correspond to the model velocity curve 
while the crosses correspond to the data. The agreement is very good along the line-of- nodes. The mismatch in data and model for radii 
beyond r~ ±1'.'7 is likely due to the fact that the inclination angle of the modeled gas disk is not well represented in the outermost parts. 
This has no impact on the derived black hole mass 

bits, we fit the geometry of the ga s disk with a set 
of til ted rings, using kinemetry (|Krainovic et all 
I2006D . The mean P.A. is 155° with twists of ±14°. 
The median inclination angle is ~ 45° ± 12°. 

• We construct a tilted-ring model for the central 
H2 gas motions, in the combined potential of the 
stars and the black hole. We account for the veloc- 
ity dispersion via a pressure term in the isotropic, 
axisymmetric Jeans' equation. The geometry of 
the disk is fixed via the kinemetry fit to the veloc- 
ity field, but we allow the overall inclination angle 
to vary. Our best-fit model has a mean inclina- 
tion angle of 34° ± 4°, and a black hole mass of 
M BH = (4.5±i £) x 10 7 M Q . We find that the incli- 
nation angle, derived through the kinemetry fit is 
not able to reproduce the overall velocity field. We 
conclude that the optimal choice would be to si- 
multaneously fit the disk geometry and the black 
hole mass. 

As our dynamical model is by construction sym- 
metric about the center, it fails to reproduce the 
asymmetry in the H2 data. When symmetrising 
the velocity field, the agreement of data and model 
is excellent. 

• Our black hole mass measurement is somewhat 
smaller but still consistent with previous gas mea- 
surements using ionized gas species. It is substan- 
tially smal l er tha n the stellar dynamical model of 
ISilge et all (|2005f ). However, it is in perfect agree- 
ment with the stellar dynamical measurement of 
Cappellari et al. (in prep.) who extract the stellar 




Fig. 17. — Black hole mass (Mbh) vs. stellar vel ocity dispersion 
along with the best fit correlation as compiled by Trcmainc et al. 
(2002). Mass measurements based on stellar kinematics are de- 
noted by circles, on gas kinematics by triangles, and on maser 
kinematics by asterisks. Filled diam onds denote mass measure- 
ments of Cen A's black hole bv lMarconi et alj l|2001fl (red), HN+06 
(blue), and this work (green). The dashed lines show the lcr con- 
fidence levels on the best-fit correlation. Our high resolution two- 
dimensional kinematic maps and dynamical modeling bring Cen A 
in full agreement with the Mbh _<t relation 



the major part of molecular hydrogen is well set- 
tled in the total gravitational potential + ^bh, 
and is therefore a good tracer of the central black 
hole mass. 

Under the assumption of gas moving on circular or- 
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kinematics from the same SINFONI data set. 

• Our black hole mass determination brings 
Cen A in agreement with the Mrh — o relatio n 
(jFerrarese fc Merrittll2000l : iGebhardt et al.ll2000ft . 

ACKNOWLEDGMENTS 

We thank the Paranal Observatory Team for the sup- 
port during the observations. We are grateful to Davor 



Krajnovic for making the Kinemetry software publicly 
available. We thank Martin Hardcastle for kindly pro- 
viding unpublished VLA data. This work has been made 
possible through support by the Netherlands Research 
School for Astronomy NOVA. NN acknowledges support 
from the Christiane-Niisslein-Volhard Foundation. MC 
acknowledges support from a PPARC Advanced Fellow- 
ship (PP/D005574/1). 



REFERENCES 



Barth, A. J., Sarzi, M., Rix, H.-W., Ho, L. C, Filipponko, A. V., 

& Sargent, W. L. W. 2001, ApJ, 555, 685 
Begeman, K. G. 1987, Ph.D. Thesis, Groningen University 
Bertola, F., Cappellari, M., Funes, J. G., Corsini, E. M., Pizzclla, 

A., & Vega Beltran, J. C. 1998, ApJ, 509, L93 
Binney, J., & Tremaine, S. 1987, Princeton, NJ, Princeton 

University Press, 1987, 747 p. 
Bonnet, H., et al. 2003, Proc. SPIE, 4839, 329 
Bonnet, H. et al. 2004, The ESO Messenger 117, 17 
Cappellari, M., Verolme, E. K., van der Marel, R. P., Kleijn, 

G. A. V., Illingworth, G. D., Franx, M., Carollo, C. M., & de 

Zeeuw, P. T. 2002b, ApJ, 578, 787 
Cappellari, M. 2002a, MNRAS, 333, 400 
Cappellari, M. & Emsellem, E. 2004, PASP, 116, 138 
Carilli, C. L., Perley, R. A., & Dreher, J. H. 1988, ApJ, 334, L73 
Clarke, D. A., Burns, J. O., & Norman, M. L. 1992, ApJ, 395, 444 
Eisenhauer, F., et al. 2003a, Proc. SPIE, 4841, 1548 
Eisenhauer, F., et al. 2003b, The Messenger, 113, 17 
Emsellem, E., Monnet, G., & Bacon, R. 1994, A&A, 285, 723 
Evans, D. A., Kraft, R. P., Worrall, D. M., Hardcastle, M. J., Jones, 

C, Forman, W. R., & Murray, S. S. 2004, ApJ, 612, 786 
Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 
Ferrarese, L., & Ford, H. 2005, Space Science Reviews, 116, 523 
Ferrarese, L., Mould, J. R., Stetson, P. B., Tonry, J. L., Blakeslee, 

J. P., & Ajhar, E. A. 2007, ApJ, 654, 186 
Gebhardt, K., et al. 2000, ApJ, 539, L13 
Gerhard, O. E. 1993, MNRAS, 265, 213 
Grandi, P., et al. 2003, ApJ, 593, 160 
Haring, N., & Rix, H.-W. 2004, ApJ, 604, L89 

N. Haring-Ncumayer, M. Cappellari, H.-W. Rix, M. Hartung, M. A. 

Prieto, K. Meisenheimer, R. Lenzen, 2006, ApJ, 643, 226 
Hardcastle, M. J., Worrall, D. M., Kraft, R. P., Forman, W. R., 

Jones, C, & Murray, S. S. 2003, ApJ, 593, 169 
Israel, F. P., van Dishoeck, E. F., Baas, F., Koornneef, J., Black, 

J. H., & de Graauw, T. 1990, A&A, 227, 342 
Israel, F. P. 1998, A&A Rev., 8, 237 
Jones, D. L., et al. 1996, ApJ, 466, L63 
Kormendy, J., & Richstone, D. 1995, ARA&A, 33, 581 
Krajnovic, D., Cappellari, M., de Zeeuw, P. T., & Copin, Y. 2006, 

MNRAS, 366, 787 
Krajnovic, D., Sharp, R., & Thatte, N. 2007, MNRAS, 374, 385 
Krause, M. 2005, A&A, 431, 45 

Macchetto, F., Marconi, A., Axon, D. J., Capetti, A., Sparks, W., 
& Crane, P. 1997, ApJ, 489, 579 



Marconi, A., Capetti, A., Axon, D. J., Koekemoer, A., Macchetto, 

D., & Schreier, E. J. 2001, ApJ, 549, 915 
Marconi, A., & Hunt, L. K. 2003, ApJ, 589, L21 
Marconi, A., Pastorini, G., Pacini, F., Axon, D. J., Capetti, A., 

Macchetto, D., Koekemoer, A. M., & Schreier, E. J. 2006, A&A, 

448, 921 

Nicholson, R. A., Bland-Hawthorn, J., & Taylor, K. 1992, ApJ, 
387, 503 

Peng, E. W., Ford, H. C, Freeman, K. C, & White, R. L. 2002, 

AJ, 124, 3144 
Pringle, J. E. 1996, MNRAS, 281, 357 

Quillen, A. C, de Zeeuw, P. T., Phinney, E. S., & Phillips, T. G. 

1992, ApJ, 391, 121 
Quillen, A. C, Graham, J. R., & Frogel, J. A. 1993, ApJ, 412, 550 
Quillen, A. C, Brookes, M. H, Keene, J., Stern, D., Lawrence, 

C. R., & Werner, M. W. 2006, ApJ, 645, 1092 
Quillen, A. C, Bland-Hawthorn, J., Brookes, M. H, Werner, 

M. W., Smith, J. D., Stern, D., Keene, J., & Lawrence, C. R. 

2006, ApJ, 641, L29 
Reipurth, B., Heathcote, S., Morse, J., Hartigan, P., & Bally, J. 

2002, AJ, 123, 362 
Rejkuba, M. 2004, A&A, 413, 903 
Schreier, E. J., et al. 1998, ApJ, 499, L143 

Shapiro, K. L., Cappellari, M., de Zeeuw, T., McDermid, R. M., 

Gebhardt, K., van den Bosch, R. C. E., & Statler, T. S. 2006, 

MNRAS, 370, 559 
Schoenmakers, R. H. M., Franx, M., & de Zeeuw, P. T. 1997, 

MNRAS, 292, 349 
Silge, J. D., Gebhardt, K., Bergmann, M., & Richstone, D. 2005, 

AJ, 130, 406 

Taylor, D., Dyson, J. E., & Axon, D. J. 1992, MNRAS, 255, 351 
Tingay, S. J., et al. 1998, AJ, 115, 960 

Tonry, J. L., Dressier, A., Blakeslee, J. P., Ajhar, E. A., Fletcher, 
A. B., Luppino, G. A., Metzger, M. R., & Moore, C. B. 2001, 
ApJ, 546, 681 
Tremaine, S., et al. 2002, ApJ, 574, 740 
van der Marel, R. P., & Franx, M. 1993, ApJ, 407, 525 
van der Marel, R. P., & van den Bosch, F. C. 1998, AJ, 116, 2220 
Verdoes Kleijn, G. A., van der Marel, R. P., de Zeeuw, P. T., Noel- 
Storr, J., & Baum, S. A. 2002, AJ, 124, 2524 



